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ABSTRACT 


Binary  liquid- metal  solubility  equilibria  of  the  type  A  (in 
solid  or  liquid  phase)  (in  liquid  b)  are  analyzed  in  terms  of 

the  fundamental  thermodynamic  solution  parameters  ah“,  and 
ASa^)>  the  heat  and  entropy,  respectively,  of  solution  of  the  liquid 
solute  at  infinite  dilution  in  the  liquid  solvent.  A  general  method 
for  the  treatment  of  equilibrium  solubility  data  is  developed  and 
is  termed  the  X'  -function  method,  in  view  of  its  similarity  to  the 
standard  x -function  method  for  the  treatment  of  vapor-pressure 
data.  By  the  assumption  of  subregular  behavior  of  the  liquid  metal 
solution,  the  X'  -function  method  is  applicable  to  moderately  con¬ 
centrated  solutions.  Consideration  is  givento  details  of  the  appli¬ 
cation  of  the  method  to  the  three  general  states  of  the  solute 
phase:  (a)  pure  solute  A,  (b)  terminal  solution  of  b  in  a,  and  (c) 
intermetallic  compound  A^i^ .  The  X' -function  method  is  applied 
to  five  binary  systems  for  which  data  are  available  in  the  litera¬ 
ture:  Mg-Fe,  Pb-Ni,  Bi-Be,  Pb-Zn,  andPb-Cu,  The  results  are  in 
close  agreement  with  the  results  published  in  the  literature. 
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A  THERMODYNAMIC  ANALYSIS  OF  SOLUBILITY  IN  LIQUID  METALS 


INTRODUCTION 

A  major  problem  in  the  efficient  utilization  of  nuclear  energy  is  the  extraction  of 
heat  from  a  reactor  operating  at  a  high  temperature.  The  properties  of  many  liquid 
metals  recommend  their  use  as  heat  transfer  fluids,  but  a  serious  obstacle  to  their 
application  lies  in  the  phenomenon  of  mass  transfer  (1).  In  general,  a  liquid  metal  cir¬ 
culating  in  a  heat  exchanger  tends  to  dissolve  the  container  metal  in  the  hot  zone  and  to 
transfer  this  metal  to  a  cooler  zone  where  it  may  deposit  and  plug  the  conduit.  In  recent 
studies  (2,3)  attention  has  been  directed  to  the  fundamental  thermodynamic  condition  for 
mass  transfer,  namely,  that  the  solubility  of  the  solid  container  metal  in  the  liquid  metal 
must  vary  with  temperature.  Although  kinetic  factors,  such  as  diffusivities  and  rates  of 
dissolution,  nucleation,  and  growth,  play  important  roles  in  determining  the  rate  of  mass 
transfer,  it  is  recognized  that  the  basic  driving  force  for  the  process  lies  in  the  varia¬ 
tion  of  solubility  with  temperature. 

In  previous  studies  of  solubility  data  for  liquid  metal  systems  (2,3)  it  was  observed 
that  the  temperature  coefficients  of  solubility  of  solid  metals  in  dilute  liquid  metal  solu¬ 
tions  correlate  well  with  the  disparities  in  sizes  of  the  solvent  and  solute  atoms.  This 
correlation  is  shown  in  Fig.  1.  The  temperature  coefficients  are  expressed  in  the  form 
-d(R  In  )/d(  1/T) ,  where  nA(^)  is  the  atom  fraction  of  the  solute  A  in  the  saturated 
liquid  metal  B,  R  is  the  gas  constant,  and  T  is  the  Kelvin  temperature.  In  a  recent  com¬ 
munication  (3),  the  thermodynamic  basis  for  this  form  of  the  temperature  coefficient  was 
briefly  outlined,  and  some  structural  implications  of  the  correlation  with  the  atomic  size 
factor  were  discussed. 

The  purpose  of  the  present  report  is  to  examine  more  rigorously  the  solubility  rela¬ 
tions  for  binary  metal  systems  and  to  present  a  method  for  the  treatment  of  solubility 
results  in  terms  of  the  basic  thermodynamic  parameters  of  the  systems. 

The  term  solubility  as  used  herein  refers  to  the  equilibrium 

A  (in  an  A-rich  phase,  usually  solidJiF^A  (in  liquid  B)  (1) 

where  the  extent  of  solution  of  A  in  liquid  Bis  determined  by  the  equivalence  of  the 
chemical  potentials  of  each  component  in  the  two  phases.  Two  important  thermodynamic 
properties  of  the  liquid  solution,  which  can  be  evaluated  from  measurements  of  such 
equilibria,  are  the  energy  and  entropy  changes  accompanying  the  transfer  of  a  solute 
atom  from  pure  liquid  solute  to^ure  liquid  solvent.  These  quantities  are  expressed  here 
by  the  partial  molal  enthalpy  AhA(^)  and  the  partial  molal  excess  entropy  of 
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ATOMIC  SIZE  FACTOR.  RATIO  OF  RADfUS  OF  SOLVENT  ATOM  TO  RADIUS  OF  SOLUTE  ATOM 


Fig.  1  -  Temperature  coefficients  of  solubility  in  binary  metallic  systems  as  a  function 
of  size  factor.  Double  symbols  (®  or  [oj  )  refer  to  systems  in  which  the  solute-rich  phase 
is  essentially  pure  solute. 
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component  A  at  infinite  dilution  in  liquid  b.*  The  method  developed  below  for  the  treat¬ 
ment  of  solubility  data  in  terms  of  ah”  and  as^>“  will  be  termed  the  S -function 
method  owing  to  its  similarity  to  one  of  the  standard  methods  for  the  treatment  of  vapor 
pressure  data  (4).  The  similarity  results  from  the  fact  that  in  both  cases  the  process  is 
one  of  the  transfer  of  an  atom  from  a  condensed  phase  to  an  attenuated  phase,  viz.,  a  gas 
in  the  case  of  vaporization  and  a  dilute  solution  in  the  case  of  solubility.  The  k- function 
method  is  directly  applicable  to  dilute  solutions  and,  by  means  of  approximations  based 
on  the  subregular  solution  model  (5) ,  t  may  be  extended  to  treat  moderately  concentrated 
solutions. 


From  the  viewpoint  of  a  further  understanding  of  alloying  behavior  (in  solid  as  well 
as  liquid  alloys),  measurements  of  solubility  are  valuable  for  two  reasons.  First,  values 
of  ah"  and  as*8’"  can  be  determined  for  systems  whose  ranges  of  composition  are  so 
low  that  ordinary  thermodynamic  methods  for  their  measurement  are  impractical.  Sec¬ 
ond,  systems  of  low  solubility  are  of  considerable  theoretical  interest  in  studying  the 
effect  of  the  atomic  size  factor  since  these  systems  commonly  have  size  effects  so  large 
as  to  overshadow  other  possible  effects  (3).  Solubility  data  for  liquid  metals  appear 
abundantly  in  the  literature,  but  in  only  a  relatively  few  cases  (e.g.,  Kleppa  (7,8))  have 


the  data  been  treated  to  evaluate  AH 


A(-t) 


and 


ASXV  and  their  corresponding  uncertainties. 

A(-t) 


DEVELOPMENT  OF  THE  GENERAL  SOLUBILITY  RELATION 

The  approach  is  to  relate  the  properties  of  the  solute  A  dissolved  in  liquid  B  to  the 
properties  of  pure  A  in  the  liquid  state.  In  most  cases  this  reference  state  of  the  pure 
liquid  solute  is  hypothetical,  at  the  temperatures  of  interest,  and  the  properties  of  the 
pure  liquid  must  be  defined  by  a  reasonable  and  continuous  extrapolation  of  the  proper¬ 
ties  above  the  melting  point. 

The  general  solubility  relation  is  developed  from  consideration  of  the  reaction 
scheme  shown  in  Fig.  2  for  component  A  of  the  system  A-B.  The  symbols  A(s)  and  A(-ti 
refer  to  the  pure  component  A  in  the  solid  and  liquid  states,  respectively.  nA(s)  and 
NA( are  the  concentrations  of  component  A  in  the  two  mutually  saturated  phases.  The 
various  terms  in  af  represent  the  tree-energy  changes  involved  in  the  transfer  of  one 
mole  of  component  A  from  one  phase  or  state  to  another,  e.g.,  af"  is  the  free  energy  of 
melting  pure  A  at  the  temperature  of  the  solubility  equilibrium. 


*The  partial  molal  excess  quantities  represent  deviations  from  the  properties  of  ideal 
solutions;  thus, 


as*{,  =  aV>  -  Asi?ir  '  av,  •  *  Vt, •  (2) 

and 

=  A?V)  -  =  aFA(*)  -  RT  -tn  ,  (3) 

_ xs 

where  AFa^^  is  the  partial  molal  excess  free  energy.  From  these  definitions  it  should 
be  noted  that 


“i’*,  =  %*)  •  (4) 

t  The  subregular  solution  model,  whose  name  was  proposed  by  Hardy  (6),  differs  from  the 
well-known  regular  solution  model  in  that  several  additional  parameters  are  available 
to  provide  for  the  deviation  of  real  solutions  from  regular  behavior. 
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LIQUID  STATES 


SOLID  STATES 


A  (AT  N 


'A(i) 


NA(j)  naU) 

-Na(ATOM  FRACTION  OF  COMPONENT  A) 


IN  LIQUID  B) 


Fig.  2  -  Reaction  scheme  for 
component  A  of  a  system  A-B, 
illustrating  the  states  and 
free-energy  changes  employed 
in  developing  the  general  sol¬ 
ubility  relation 


The  solid  and  liquid  phases  are  at  equilibrium  with  each  other;  therefore  AFj  =  o  , 
and  (see  Fig.  2) 


AF 


Apt) 


+  AF,  -  AF 


A(s) 


(5) 


Since  the  thermodynamic  properties  of  the  liquid  binary  solution  are  determined  by  AfA(,^ 
and  its  variation  with  temperature,  the  procedure  is  to  evaluate  each  free-energy  change 
in  Eq.  (5)  in  terms  of  NA(^  .  NA(s)>  A(-t)  >  asa(-E)>  and  the  other  relevant  thermody¬ 
namic  properties  of  the  system  A-B  . 

The  free  energy  of  melting  of  pure  A,  AFa,  can  be  obtained  by  integration  of 


-  as"1 

dT 


(6) 


where  asa  is  the  entropy  of  melting  of  pure  A  at  temperature  t.  asa  iB  given  by 


AS. 


AS 


A,T' 


u 


feu 


p  •  ^  I 

T  ~  dTj 


(7) 


where  ta  is  the  melting  point  of  pure  A  and  ac™  A  is  the  difference  in  heat  capacities  of 
liquid  and  solid  A.  Integrating  Eq.  (6)  and  making  use  of  the  fact  that  af"1  „  =  0 


A.T. 


,  mm  m 

Afa  =  ASA  T-  (Ta  -  T) 


-r  r 


AC. 


P,  A 


dT  dT  . 


(8) 


However,  if  component  A  undergoes  transitions  in  the  solid  state  or  if  AC™  A  cannot 
be  represented  by  a  simple  analytic  function,  the  integral  in  Eq.  (8)  becomes  cumbersome 
to  evaluate  and  it  is  more  satisfactory  to  make  use  of  the  free-energy  functions  to  eval¬ 
uate  Afa.  The  free-energy  function,  fef ,  is  defined  (9)  by 


fef 


,  (%^) .  St 


(9) 


where  T  is  the  Kelvin  temperature,  st  refers  to  the  standard  state  at  a  specified  stand¬ 
ard  temperature  (usually  298.15°K),  and  ST  is  the  absolute  entropy  at  T°K .  The  free- 
energy  functions  are  readily  constructed  for  the  solid  and  liquid  over  their  respective 
temperature  ranges  of  stability  by  evaluating  the  terms  on  the  right-hand  side  of  Eq.  (9) 
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from  standard  tables  (9  to  14).  An  important  practical  characteristic  of  the  free-energy 
function  is  that  it  varies  relatively  slowly  with  temperature,  so  that  a  plot  of  the  free- 
energy  function  as  a  function  of  temperature  for  the  condensed  phases  of  a  material  can 
be  readily  constructed  from  a  relatively  small  number  of  points. 

The  construction  of  the  free-energy  function  for  the  liquid  at  temperatures  below  the 
melting  point  is  usually  accomplished  by  assuming  that  a  simple  heat  capacity  expres¬ 
sion,  e.g., 


CP. A(t) 


a  +  bT  , 


(10) 


is  representative  of  the  liquid  properties  below,  as  well  as  above,  the  melting  point.  If 
Eq.  (10)  is  used,  the  free-energy  function  of  pure  liquid  A  is  then  given  by 


+  a 


where 


b 

+  2 


(11  .  -  H 


St)A(.f) 


and  s 


tT.a(*) 


(11) 


are  the  enthalpy  (relative  to  the  standard  temperature  st)  and  the  absolute  entropy  of  the 
pure  liquid  A  at  its  melting  point.  Since  the  free-energy  functions  for  solid  and  liquid  A 
both  relate  to  the  same  standard  state  of  pure  A  at  the  standard  temperature  st,  the  free 
energy  of  melting  is  given  by 


(12) 


The  term  afA(  s)  of  Eq.  (5)  provides  for  the  fact  that  the  activity  of  A  in  the  phase  in 
equilibrium  with  the  liquid  solution  may  be  lower  than  that  of  pure  A .  Formally,  this 
term  is  given  by 

^A(,s)  "  RT  tn  aA(s)  =  RT  -fn  {>'A(  s)  NA(  s))  (13) 

where  aA(s),  rA(s)>  and  na(S>  are  the  activity, _activity  coefficient,  and  atom- fraction  of 
component  A  in  the  solid  phase,  respectively.  AFA(.  may  be  evaluated  exactly  if  activi¬ 
ties  and  phase  boundaries  have  been  well  determined  for  the  solid  phase;  if  these  data 
are  not  available,  resort  must  be  made  to  the  approximations  to  be  treated  later. 


The  term  af A(^  of  Eq.  (5)  may  be  written  as 

A*A(-t)  =  AHA(,}  -  T  AS^  ; 


(14) 


AF 


A(-t) 


ah 


AC-t ) 


+  RT  In  N 


A(-t) 


-xs 

ASA(-£) 


(15) 


ASA(f.)  represents  the  deviation  from  the  ideal  entropy  of  mixing  calculated  on  a  substi¬ 
tutional  basis  and  may  consist  of  thermal  (i.e.,  vibrational)  as  well  as  configurational 
contributions.  aha^  represents  the  energy  difference  between  A  atoms  in  the  liquid  B 
solution  and  A  atoms  entirely  surrounded  by  A -neighbors,  as  in  pure  liquid  A. 
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The  general  solubility  relation  is  obtained  by  combining  Eqs.  (5),  (13),  and  (15);  thus, 

+  RT  In  NA(  t)  -  T  AS**^  +  AFa  -  RT  -Cn  aft(  s)  .  (16) 

Regrouping  and  dividing  by  t, 


afa 

R  -fn  NaU)  +  -y  -  R  aA( s ) 


AH 


MAl  +  Aoxs 

T  AbA  U) 


(17) 


In  this  form  the  solubility  relation  is  suitable  for  application  to  the  several  types  of 
solubility  equilibria.  For  convenience,  let  the  left-hand  side  of  Eq.  (17)  be  designated  as 
the  ^-function,  i.e., 


X  =  R  NA(^} 


AF.  -  AF 


A(s) 


=  R  -tn 


\aA(s)  j 


(18) 


In  general,  and  asa(s^  are  functions  of  temperature  and  composition,  but  for 

very  dilute  solutions,  e.g.,  nA(  }  <  0.005,  the  composition  dependence  may  be  neglected. 
Thus,  if  the  values  of  X  calculated  from  experimental  results  by  Eq.  (18)  are  plotted 
against  the  reciprocal  of  the  Kelvin  temperature,  the  slope  of  the  tangent  to  the  resulting 
curve  is  -AH^(^)  and  the  intercept  of  the  tangent  at  1/T  -  0  is  asa*|".  Since  AH^  may 
be  a  function  of  temperature,  the  x  values  will  generally  define  a  gentle  curve  rather 
than  a  straight  line.  However,  in  actual  cases,  strong  deviations  from  the  Kopp- Neumann 
rule  for  the  additivity  of  heat  capacities  are  not  expected  for  liquid  alloys,  and,  in  view  of 
the  scatter  to  be  expected  in  experimental  results  for  the  solubilities,  a  straight  line  usu¬ 
ally  offers  the  most  satisfactory  fit. 


In  many  binary  metal  systems  the  eutectic  occurs  at  solute  concentrations  greater 
than  nA(,£)  =  0.005,  and  it  is  most  desirable  to  extend  the  X -function  method  to  treat 
these  more  concentrated  solutions.  Lumsden's  extension  of  the  regular  solution  model 
(15),  termed  the  subregular  solution  model  (5,6),  has  been  found  to  give  a  good  repre¬ 
sentation  of  the  thermodynamic  properties  of  many  binary  liquid  systems  (10)  and  of.ers 
a  convenient  method  of  introducing  a  first  approximation  for  the  concentration  depend¬ 
ence  of  AHAa)  and  asa^j.  Thus,  let 


AHA (•£)  ~  AHA(-t)  ~  AHA(£>  ~ 


and 


ASA(^)  =  ASA(-t)  NB(^)  ASA(i)  (1  ~ 


Substituting  in  Eq.  (17), 


x  = 


“  na^)> 


_ 03 

_  -  _  AHA(t)  .  a-bxs'“ 

2  "  T  aA(-£)  • 


(19) 

(20) 


(21) 


Values  of  X'  calculated  from  experimental  results  for  more  concentrated  solutions  may 
now  be  plotted  against  l/T,  and  ahA(,^,  and  asa^?  are  obtained  from  the  slope  of  the 
resulting  curve  and  its  intercept  at  l/T  =  o.  Although  the  precision  of  this  treatment  for 
moderately  concentrated  solutions  depends  on  how  closely  a  given  liquid  alloy  conforms 
to  the  subregular  behavior  postulated  by  Eqs.  (19)  and  (20),  the  method  appears  to  be  fre¬ 
quently  applicable  to  concentrations  as  high  as  nA(^)  =o,i,  which  includes  most  systems 
of  interest  to  current  studies. 
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APPLICATION  OF  THE  S- FUNCTION  METHOD  TO 
THREE  GENERAL  CASES 

In  calculating  the  X  -function  by  Eq.  (18),  data  are  frequently  not  available  to  eval¬ 
uate  the  term  afa  =  RT  -tnaA(sj ,  and  approximations  are  then  necessary  in  applying 
the  I -function  method.  Thus,  it  is  convenient  to  distinguish  three  types  of  solubility 
equilibria  according  to  the  nature  of  the  phase  in  equilibrium  with  the  liquid  solution. 

The  cases  illustrated  in  Fig.  3  are:  I,  the  A -rich  phase  has  negligible  solubility  for  B ; 
H,  the  A -rich  phase  is  a  terminal  solid  or  liquid  solution  showing  appreciable  solubility 
for  b;  in,  the  A-rich  phase  is  an  intermetallic  compound. 

It  should  be  noted  that  a  distinction  between  solid  and  liquid  A  -rich  phases  is  unnec¬ 
essary.  If  the  A-rich  phase  is  a  liquid,  then  Afa  -  0  in  the  calculation  of  the  2-function. 
If  the  free-energy-function  method  is  employed  to  evaluate  afa,  the  phase  transitions  in 
A  are  properly  accounted  for  in  the  calculation  of  afa,  and  all  solubility  data  may  be 
treated  in  one  plot,  irrespective  of  the  number  of  transformations  or  reactions  of  the  A- 
rich  phase  over  the  temperature  range  of  the  solubility  data. 


Case  I 

If  the  A-rich  phase  shows  a  maximum  solubility  for  B  of  two  atomic  percent,  or 
less,  the  contribution  of  afa,  s .  to  1  will  generally  be  less  than  one  percent.  Except  for 
the  most  precise  solubility  data,  afA(  s)  may  then  be  neglected,  and  the  2-function 
becomes 


2  -  R  tn  NA(-£5 


(22) 


Case  H 

In  cases  of  appreciable  solubility  of  component  B  in  the  A  -rich  terminal  solution  it 
is  frequently  found  that  the  phase  boundaries  in  equilibrium  with  the  liquid  solution  are 


Case  I  Case  II  Case  HI 

Fig.  3  -  Three  cases  of  liquid  metal  solubility  equilibria  distinguished 
according  to  the  type  of  approximation  employed  in  evaluating  the 
activity  of  the  solute  A  in  the  A-rich  phase 
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at  least  approximately  known  but  that  the  thermodynamic  data  necessary  to  evaluate  yh(  s) 
are  not  available.  Except  in  cases  of  unusually  high  solubility  (i.e.,  nA(  s)  >  o.i),  it  is 
reasonable  to  evaluate  af  A(  s )  by  the  assumption  of  Raoult's  law;  thus,  -yA(s)  =  l  and 


2  -  R  -tn 


(23) 


Case  m 


For  systems  in  which  the  liquid  solution  is  in  equilibrium  with  an  intermetallic 
phase,  the  thermodynamic  information  necessary  for  an  exact  evaluation  of  afa^  s)  will 
only  infrequently  be  available.  In  the  relatively  rare  cases  where  the  composition  range 
of  the  intermetallic  phase  is  very  broad,  e.g.,  the  phase  y3'  in  Case  in  of  Fig.  3,  such 
thermodynamic  data  are  indispensable  for  an  evaluation  of  the  liquid  solubility  in  terms 
of  ahA(^)  and  asa“^.  However,  in  most  cases,  the  mtermetallic  phase  is  an  ordered 
compound  of  very  limited  composition  range,  and  afA(  s)  may  be  approximated  if  the  free- 
energy  function  for  the  intermetallic  compound  is  known. 


For  the  formation  of  an  intermetallic  compound  A^  _  x  from  the  pure  components, 

x  A  +  (1 -x)  b  ^  AXB1.X ,  (24) 

A  fef  (AXB1^X)  =  fefA  -  [xfefA  +  (l-x)fefg]  .  (25) 


The  free  energy  of  formation  of  AxBi.x  is  given  by 


AFf(AxBi.x)  AH^A^.J 


+  A  fef  (AxBj.x) 


(26) 


If  the  free-energy  functions  for  A,  B,  and  are  known,  and  if  the  standard  heat  of 

formation  of  AxBl  x  has  been  measured,  then  AFf(AxBi.x)  can  be  calculated  for  any 
temperature. 

The  method  of  estimating  afA(s)  from  AFf(AxBi.x)  is  indicated  in  Fig.  4.  The  par¬ 
tial  molar  free  energy  of  component  A  in  the  intermetallic  compound  in  equilibrium  with 
the  liquid  solution  is  given  by  the  intercept  on  the  axis  na  =  l  of  the  tangent  common  to 
the  free-energy  curves  of  the  intermetallic  compound  and  the  liquid  solution.  In  cases  of 
limited  solubility,  the  free-energy  curve  of  the  liquid  solution  must  turn  sharply  upward 
within  the  first  few  atom  percent  of  solute  concentration.  Thus,  only  a  second-order 
error  is  introduced  by  placing  the  tangent  line  through  the  origin  rather  than  tangent  to 
the  as-yet-unknown  curve  for  the  liquid  solution.  The  magnitude  of  the  error  will  depend 
on  the  shape  and  position  of  both  free-energy  curves;  however,  it  is  clear  that  the  error 
is  less  for  those  systems  in  which  (a)  the  free-energy  curve  for  the  intermetallic  com¬ 
pound  is  sharper  and  deeper,  (b)  the  liquid  solubility  is  more  limited,  and  (c)  the  propor¬ 
tion  of  component  A  in  the  intermetallic  compound  is  greater. 

Using  the  approximation  for  the  tangent  line, 


The  2-function  then  becomes 


AF 


A(s) 


AFf(AxBt.x) 

X 


(27) 


AFa  AFf(A^.,) 
A(-J)  T  *T 


2  “  R  -in  N 


(28) 
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Fig.  4  -  Free-energy  diagramfor  a  system 
A-B  containing  an  intermetallic  compound 
A,Bl.  xin  equilibrium  with  a  B~rich  liquid 


Note  that  any  phase  transformations  in  the  intermetallic  compound  are  properly  dealt 
with  through  the  use  of  the  free- energy- function  method  for  the  calculation  of  AFf(AxBi.*). 


EVALUATION  OF  THE  CONSTANTS  IN  LOG  s  =  C/T  +  D 
A  common  expression  employed  to  fit  solubility  results  is 

log  s  -  +  D  (29) 


where  the  solubility  s  may  be  expressed  in  a  variety  of  units  and  c  and  D  are  assumed 
to  be  constants,  c  corresponds  to  the  temperature  coefficient  of  solubility  employed  in 
earlier  correlation  studies  (2,3).  It  is  of  considerable  interest  to  evaluate  C  and  D  in 
terms  of  the  thermodynamic  solution  parameters  and  AS  appearing  in  the  £- 

function  treatment,  cf.  Eq.  (21).  For  this  purpose,  Eq.  (29)  is  rewritten  as 


log  s 


'tn<kNAft))  C 
■in  10  T 


(30) 


where  k  is  the  ratio  of  the  solubility  s ,  in  whatever  units  have  been  used,  to  the  atom- 
fraction  Na,p-;  i,p. , 

'H-w 


Thus. 


■in  N 


A(-t) 


k  =  -sj-A— 

NA<4) 
2.303C 


+  2.  303D  -  In  k  . 


(31) 

(32) 
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To  establish  a  parallel  relation  from  the  I -function  treatment,  afa  and  afa(S)  are 
first  written  explicitly  as  ' 


mm  m 

AF.  =  AH.  -  T  AS.  , 


(33) 


^AC)  =  AHA(S)  +  ™  NA(s)  "  T  ASA(s)  • 


(34) 


Equation  (21)  may  now  be  written  as 
X  =  R  In  N 


AHA  m  AHA, s .  »  — xs 

A(-£)  +  ~ f  ASA  ~  T  R  NA(s)  +  ASA(s) 


5a  a) 


(35) 


Rearranging, 


in  NA<(f)  = 


(1  ~NA(i))2  AH^)  *  AHa  -  AHA(s.) 
RT 


(l-NAa))  AS 


A  U) 


+  AS. 


— xs 
ASA( s ) 


+  in  N 


A(  s) 


(36) 


Restricting  attention  to  dilute  solutions,  the  terms  c  and  d  may  be  identified  by  com¬ 
parison  of  Eqs.  (32)  and  (36);  thus, 


AH 


C  =  - 


-VU  III  - 

[A(i)  +  AHA  "  ^AC  s) 
2. 303  R 


AS* 


D  = 


A(i) 


.  m  -xs 
+  ASa-ASAis) 


2.303  R 


+  log  (kNA(B)) 


(37) 

(38) 


These  expressions  for  c  and  D  are  rigorously  correct  only  when  the  various  terms 
in  AH,  AS,  and  NA(S.  are  independent  of  temperature.  In  many  cases  these  terms  vary 
relatively  slowly  with  temperature  so  that,  for  dilute  solutions,  it  is  often^possible  to 
apply  Eqs.  (33)  and  (38)  to  obtain  approximate  values  of  ah”(,^  and  ASa!^  from  a  plot  of 
log  s  versus  l/T.  However,  such  plots  should  show  abrupt  changes  in  the  slope  and 
intercept  at  temperatures  of  transition  or  with  other  changes  in  the  nature  of  the  A-rich 
phase.  These  changes  in  slope  and  intercept  are  proportional,  respectively,  to  the 
enthalpy  and  entropy  changes  of  component  A  in  the  A-rich  phase. 


THE  THERMODYNAMIC  CONDITIONS  FOR  LOW  SOLUBILITY 

In  view  of  toe  unusual  behavior  of  several  systems  for  which  the  solvent- to- solute 
radius  ratio  is  large  (2,3),  it  is  of  interest  to  examine  Eq.  (36)  to  see  what  values  of 
AHA(^j  and  As*(s^'  correspond  to  very  low  solubilities.  Ignoring  the  correction  term  for 
solubility  in  the  A-rich  phase  (unimportant  in  the  systems  for  which  the  unusual  behavior 
has  been  observed),  Eq.  (36)  may  be  expressed  for  dilute  solutions  as 

—co  m  — xs,  w  m 

*  „  AHA<4>  +  ^A  ASa(4)  +  ASA 

U  "»<*>  - - RT -  * - R - 


(38) 
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A  general  rule  for  the  estimation  of  the  entropy  of  fusion  of  metals  is  given  by 
Kubaschewski  and  Evans  (16)  as 


Thus, 


2.2  cal 


°K 


- 1 


- 1 

gm-atom 


m  m 

AH  m  =  2.2  T«  cal  gm-atom 


.  l 


Applying  these  approximate  relations  to  Eq.  (39), 


(40) 


(41) 


■in  N 


—00  _xs,® 

AHA(l)  asa(1) 
Hi)  '  ~  RT  R 


2.2 

R 


(42) 


The  last  term  represents  the  solubility  in  a  system  in  which  the  liquid  solution  displays 
ideal  behavior.  Thus,  the  temperature  coefficient  of  solubility  for  a  system  obeying 
Raoult's  law  and  whose  solute  follows  the  rule  expressed  by  Eq.  (37)  is  given  by 


d  (R  -in  NA(4)> 
d(l/T> 


(43) 


For  liquid  alloys  in  which  the  solute  atoms  occupy  substitutional  sites  in  the  liquid 
structure,  the  entropy  of  formation  is  expected  to  be  near  that  calculated  for  ideal  solu- 
tions.OT  In  such  cases,  which  appear  to  comprise  the  majority  of  liquid  metal  solutions, 
asa^”  is  small  (10^16,17).  For  these  systems  Eq.  (42)  shows  that  very  low  solubilities 
obtain  only  when  ahA(.p)  is  large  and  positive.  Thus  a  very  low  solubility  is  ordinarily 
accompanied  by  a  high  temperature  coefficient.  For  example,  if  ASA?,f“  =  0  and  _<n 
Ta=  1500°K,  a  solubility  of  nA(^)  =  10'5  at  T  =  1000°K  corresponds  to  a  value  of  ahA(/ 
near  22  kcal  gm-atom  ’ 1 . 


In  the  systems  demonstrating  unusual  behavior,  e.g.,  Fe  or  Ni  in  liquid  Na  (2),  very 
low  solubilities  are  found,  together  with  low  temperature  coefficients.  In  such  cases 
<£  has  a  large  negative  value.  Thus,  if  AMA(4)  =  5  kcal  gm-atom'1  andTA  =>  1500°K, 
a  solubility  of  ~  10' 5  at  1000°K  corresponds  to  a  value  of  asa^"  near  -17  cal  °K‘1 
gm-atom  • *.  Such  deviations  from  ideal  substitutional  behavior  are  too  large  to  be 
attributed  to  thermal  contributions,  which  would  involve  very  large  deviations  from  the 
Kopp-Neumann  rule.  In  terms  of  configurational  contributions,  the  most  elementary 
explanation  for  large  negative  values  of  asa*.>”  is  that  the  number  of  favorable  sites  for 
solute  atoms  in  the  solvent  matrix  is  much  smaller  than  the  total  number  of  either  sub¬ 
stitutional  or  interstitial  sites. 


EXAMPLES  OF  THE  APPLICATION  OF  THE 

2- FUNCTION  TREATMENT 

This  report  is  concluded  by  the  presentation  of  some  examples  of  the  use  of  the  £- 
function  method  for  the  treatment  of  solubility  data  reported  in  the  literature  for  several 
binary  liquid- metal  systems.  The  treatment  is  given  in  some  detail  for  the  first  system, 
Mg-Fe;  for  the  other  systems  only  the  essential  features  are  indicated.  For  brevity,  all 
numerical  values  of  thermodynamic  quantities  are  expressed  in  terms  of  the  absolute 
calorie,  the  gram-atom,  and  the  Kelvin  degree.  Thus,  the  units  of  cp  and  AS  are  cal 
“K"1  gm-atom'1,  and  those  of  ah  and  af  are  cal  gm-atom'1. 
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Magne  sium-  Iron 

The  solubility  of  solid  iron  in  liquid  magnesium  is  slight.  The  solid  phases  in 
equilibrium  with  the  liquid  solution  are  essentially  pure  Fe  (18).  Over  the  temperature 
range  for  which  solubility  data  are  available,  Fe  undergoes  two  solid-state  transitions, 
viz.,  the  a -7  transition  at  1183°K  and  the  magnetic  transition  occurring  in  the  neighbor¬ 
hood  of  1033°K  (10).  Several  good  sets  of  solubility  determinations  are  available  (19-21). 
In  each  case  the  solubility  was  determined  by  wet  chemical  analysis  of  samples  which 
had  been  held  at  controlled  temperatures  for  periods  long  enough  to  ensure  saturation  of 
the  liquid  Mg . 

In  evaluating  the  2-function  it  is  first  noted  that  the  highest  solubility  reported  is 
0.84  weight-percent  Fe,  which  corresponds  to  NFe(^)  =  0.0037.  The  correction  for  con¬ 
centrated  solutions  is  thus  small  but  not  neglected  in  the  following  treatment  of  the  data. 
Since  the  solubility  of  Mg  in  solid  Fe  is  negligible,  afFc(s)  is  zero,  and  the  2-function 
reduces  to 


2(Fe  in  Mg.-t)  =  R  in  WFe(i) 


(44) 


The  term  AFFe/T  is  evaluated  by  means  of  Eq.  (12).  The  free-energy  function  for 
Fe(Jf }  is  constructed  by  extrapolation  of  data  at  higher  temperatures  using  the  heat 
capacity  expression. 


Cp . Fe{-f)  =  9.72  +  0.43  X  10- 3T. 


(45) 


Thus,  applying  Eq.  (11), 


-f2  -  22.381  log  T  -  0.215  x  10'3T  +  59.382  , 


(46) 


where  st  is  298.15°K.  Values  of  (FT-Hst/T)Fe(s)  are  obtained  from  tables  (10);  the 
phase  transitions  in  solid  Fe  prevent  the  development  of  a  simple  analytic  expression 
similar  to  Eq.  (45).  The  values  of  afFr/t  calculated  by  Eq.  (12)  are  plotted  in  Fig.  5. 

The  changes  in  slope  at  the  transition  temperatures  are  barely  discernible  owing  to  the 
small  enthalpy  changes  accompanying  these  transitions  (e.g.,  AhJ'^  =  215  (13)  and 
AH™4c  *  410  (14)).  Values  of  afFc/T  calculated  from  Eq.  (8),  on  toe  assumption  that 
Acp  is  negligible,  are  included  in  Fig.  5  to  show  the  difference  from  the  values  calculated 
by  the  more  exact  free- energy-function  method.  At  temperatures  near  the  melting  point 
of  Fe  (1812 °K)  the  differences  are  negligible,  but  at  lower  temperatures  the  differences 
become  important,  e.g.,  0.2  cal  “K-1  gm-atonr1  at  1000°K. 

The  solubility  results  in  the  literature  have  been  reported  in  units  of  weight-percent 
of  Fe .  The  evaluation  of  2  for  each  datum  point  can  be  expedited  by  noting  that  for 
dilute  solutions 


Fe(-f)  WFe(f)  M 


Me 

Fe 


=  w. 


24.32 
fe(t)  55.  R5 


(47) 


where  WFe{^}  is  the  weight-percent  of  Fe  in  Mg(^),  and  MMg  and  MFe  are  the  atomic 
weights  of  Mg  and  Fe.  Using  Sq.  (47), 


AF 


m 

Fe 


2(Fe  in  Mg.-t  )  =  4.576  log  WFe  -  1 . 652  + 


T 


(48) 
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Fig.  5  -  The  term  AFf„/T  of  the  ^-function  for  Fe(s)  ,  The  upper  curve  is 
obtained  by  use  of  Eq.  (8)  on  the  assumption  that  ACp  is  negligible.  The  lower 
curve,  calculated  by  means  of  the  free-energy  functions  for  solid  and  liquid  Fe , 
is  more  exact  and  takes  account  of  the  transitions  in  Fe(s)  at  1183°K  and  near 
1033°K. 


Values  of  2'(Fe  in  Mg.-fc)  =  S(Fe  in  Mg,-t)/(i  -^vC(Z02  >  calculated  from  the  results  of 
Mitchell  (19).  Fahrenhorst  and  Bulian  (20),  and  Siebel  (21),  are  plotted  in  Fig.  6  as  a 
function  of  l/T.  These  data  are  replotted  on  a  reduced  scale  in  Fig.  7  to  show  the  inter¬ 
cept  at  l/T  =  o  of  the  best  straight  line  fitted  to  the  data.  The  values  of  ah”  and 
ASFe$)  Obtained  for  each  set  of  data  are  given  inTable  1.  The  2-function  plot  is  well 
suited  for  the  estimation  of  the  uncertainties  in  and  as**’T ,  which  are  also 

given  in  Table  1.  The  data  have  also  been  treated  by  plotting  log  wFe  versus  l/T .  where 
wFe  is  the  weight-percent  of  Fe  in  liquid  Mg,  and  the  results  for  ah"  and  AsF;'^  are 
also  given  in  Table  1.  In  this  case  the  two  methods  give  similar  results  owing  to  the  fact 
that  all  data  refer  to  dilute  solutions;  however,  the  2-function  method  was  found  to  be 
more  satisfactory  in  estimating  the  uncertainties  in  Ali“e(^)  and  ASpj^j. 

Weighted-mean  values  for  AHFe(^)  and  As^’^,  as  obtained  by  the  2-function 
method,  have  been  selected  as  follows: 


AHFe(^)  =  13,700  +  600, 


=  “1.5  ±  0.6  . 
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Fig.  6  -  Values  of  2'(Fe  in  Mg,£)  calculated  from 
solubility  data  reported  in  the  literature 


Lead- Nickel 

The  solubility  of  Pb  in  solid  Ni  is  believed  to  be  slight  { 1.2  atomic  percent  or 
less  (18)  and  is  therefore  neglected  in  calculating  2(Ni  in  Pb).  However,  the  solubility  of 
Ni  in  liquid  Pb  is  large,  and  measurements  have  been  reported  at  solubilities  as  high  as 
NNi  =  0.24.  Since  the  eutectic,  i.e.,  the  point  of  minimum  solubility,  lies  at  nn.  =  0.003, 
the  2-function  as  modified  for  extension  to  moderately  concentrated  solutions  is  required 
cf.  Eq.  (21). 

Two  sets  of  solubility  data  (22,23)  are  available  for  treatment,  and  values  of 
2VNi  in  Pb.-t)  are  plotted  in  Fig.  8.  The  data  of  Alden,  Stevenson,  and  Wulff  (23)  are 
given  the  greatest  weight  in  evaluating  the  thermodynamic  parameters: 

=  3050  *  600  ■ 

=  ~2-30  ±  °-6- 
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Fig.  7  -  The  2'  -Exaction  values  of  Fig.  6  plotted  on  a 
reduced  s cade  to  showthe  intercept  at  1/T  =  0.  AHFe(^  = 
13,700  ±  600;  AS*®^  =  -1.5  ±  0.6. 


Table  1 

Values  of  a5C-,j>,  and  as*3;”,  for  Mg-Fe 


AHFe(-t) 

(cal  gm-atom'1) 

AS 

(cal  °K' 1 

M  S  ,  00 

Fe(-t) 

gm-atom' 

Ref. 

Data  Treatment 
by  2-  Function 

Data  Treatment  by 

log  wFe  -  C/T  +  D 

Data  Treatment 
by  2 -Function 

Data  Treatment  by 
log  WFe  =  C/T  +  D 

13,200  ±  800 

13,400 

-1.8  ±  0.7 

-1.5 

19 

14,200  ±  700 

14,200 

-1.1  ±  0.6 

-1.2 

20 

12,600  ±  2000 

12,900 

-2.9  ±  2.0 

-2.6 

21 
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Fig.  8  -  The  5!-function  treatment  of  data  for  the  solubility  of  Ni  in  liquid  Pb. 

ARma>  ~  3050  ±  60°:  =  -2-30  ±  °-6- 


Values  of  E(Ni  in  Pb,  t) ,  not  divided  by  the  subregular  solution  factor  (l  -NNi)2,  are 
indicated  in  Fig.  8  for  the  data  of  Alden  et  al.  to  show  the  value  of  the  modified  S'- 
function  for  the  treatment  of  moderately  concentrated  solutions.  In  this  case  2'(Ni  inPb.-t) 
is  linear  as  far  as  the  solubility  data  extend,  i.e.,  to  NNi  =  0.24. 


Bismuth-  Beryllium 

Measurements  of  the  solubility  of  Be  in  liquid  Bi  have  recently  been  reported  by 
Horsley  and  Maskrey  (24).  These  workers  found  no  evidence  for  the  existence  of  inter- 
metallic  compounds;  the  values  of  2(Be  in  Bi.-C)  have,  therefore,  been  calculated  assum¬ 
ing  that  the  solid  phase  is  pure  Be.  The  values  of  AFge/T  were  constructed  from  the  data 
and  estimates  reported  by  Stull  and  Sinke  (14);  a  value  of  Cp  Qe(^)  »  7.5  was  used.  The 
values  of  the  thermodynamic  parameters  evaluated  from  the  data  shown  in  Fig.  9  are 

AHge(£)  =  22,400  ±  4000  , 

-xs.® 

&SBe(,£)  =  5.9  ±  2.5  . 


Lead- Zinc 

The  Pb-Zn  system  is  of  special  interest  since  the  metals  display  immiscibillty  in  the 
liquid  state  as  well  as  the  solid  state  (18).  Extensive  data  are  available  for  both  the 
lead-rich  and  the  zinc-rich  solubility  equilibria,  and  at  least  two  attempts  have  been 
made  to  develop  a  self-consistent  set  of  thermodynamic  relations  through  detailed  eval¬ 
uation  of  all  thermodynamic  and  phase  data  available  for  the  system  (10,25). 

The  approach  employed  here  in  applying  the  ^-function  treatment  to  the  Pb-Zn  equi¬ 
librium  data  is  based  on  the  method  of  successive  approximations  since;  in  principle,  the 
thermodynamic  results  for  the  solute  phase  are  necessary  in  order  to  evaluate  rigorously 
the  data  for  the  solvent  phase.  The  data  for  the  phase  of  more  extensive  solubility,  i.e., 
the  Pb-rich  phase,  are  evaluated  first  since  the  approximations  necessary  in  evaluating 
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Fig.  9  -  The  2-function  treatment  of  data  for  the  sol¬ 
ubility  of  Be  in  liquid  Bi .  =  22,400  ±  4000; 

C;i)  =  5^2'  5. 


_ CO  _ xs  , 00 

the  2 -functions  are  less  extensive.  On  the  basis  of  the  values  of  AH2n(Pb^)  and  AS  ^ 
derived  for  the  Pb-rich  phase,  the  2-functions  for  the  Zn-rich  phase"  are’computed,  and" 
the  corresponding  values  of  AHph(Zri  ^  and  ASpb’ Zn  ^  are  evaluated.  A  second  series 
of  evaluations  can  then  be  performed  if  the  data  warrant  it. 

2-functions  for  the  Pb-rich  phase  have  been  calculated  from  the  equilibrium  data  of 
Waring,  Anderson,  Springer,  and  Wilcox  (26),  Kleppa  (27),  and  Hass  and  Jellinek  (28),  and 
the  results  in  terms  of  2'(Zn  in  Pb.-t)  are  plotted  in  Fig.  10.  The  form  of  the  2'-function 
is 


2'(Zn  in  Pb,-£) 


1 


2 

NPb(Pb.-f  ) 


R  'tn 


(  NZn(Pb,-t) 
^ZnfZn.'t) 


(49) 


where  NZn(Zn  t,  =  1  below  TZn  =  692.6°K  and  6FZn  =  0  above  TZn.  The  data  of  Spring 
and  Romanoff’  (29)  are  inordinately  high  and  were  not  used.  Lumsden  has  evaluated  the 
eutectic  point  from  the  experimental  results  of  Hodge  and  Heyer  (30)  and  Heycock  and 
Neville  (31),  and  the  2 '-functions  are  also  shown  in  Fig.  10. 

The  values  of  the  thermodynamic  parameters  corresponding  to  the  line  shown  in 
Fig.  10  are 

^Zn(Pb^)  =  4400  ±  800 
^Zn(Pb^)  -  -0.1  ±  1.0. 
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Fig.  10  -  The  2 -function  treatment  of  data  for  the  solubility  of  Zn  in  liquid  Pb  . 
^Zn(Pb^)  =  4400  ±  800;  AS*^b  ^  =  -0.1  ±  1.0. 


The  form  of  the  2' -function  for  the  Zn-rich  phase  is 


Z'(Pb  in  Zn,  -t)  =  - - 

NZn(Zn.^) 

where  a”n(.pb  ^  is  the  subregular  solution  parameter  representing  the  deviation  from 
ideal  solution  behavior;  a£n(pb  j,)  is  given  by 


S  In 


NPb(Zn,^)\  _  NZn(Pb,l)  aZn(Pb , -t) 
i  ^Pb(Pb ,'t)/  T 


(50) 


CO  —  xs,«  _00  _XS,® 

aZn(Pb,4)  ”  AFZn(Pbl^)  “  ^ZnfPb,^)  “  T  ASZn(Pb,^)  ’  (51) 

In  view  of  the  fact  that  AS^pb  ^  is  nearly  zero,  the  calculations  can  be  simplified  by 
assuming  that  °-7n(Ph  4.  is  temperature  independent;  the  value  of  a2n(Pb  ^  was  there¬ 
fore  assumed  to  be  4450.  v  ’  ’ 


Values  of  the  2'  -functions  calculated  from  the  data  of  Waring,  Anderson,  Springer, 
and  Wilcox  (26),  Kleppa  (27),  Hess  and  Jellinek  (28),  Lumsden  (25),  Heycock  and  Neville 
(31),  and  Spring  and  Romanoff  (29)  are  plotted  in  Fig.  11.  The  thermodynamic  parame¬ 
ters  corresponding  to  the  straight  line  placed  through  the  data  are 

^Pbczn.-f)  =  12,900  ±  1300 


ASPb(Zn,-t)  "  7>  0  ±2.0. 

The  scatter  of  the  data  of  Figs.  10  and  11  indicate  that  little  improvement  is  to  be 
expected  through  another  series  of  evaluations  based  on  the  data  derived  above.  If  a  sec¬ 
ond  series  were  to  be  undertaken  it  would  be  desirable  to  employ  the  subregular  solution 
equations  (5)  to  express  the  functional  dependencies  of  the  entropies  and  heats  of  solution 
on  the  temperature  and  the  composition. 
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Fig-  11  -  The  X -function  treatment  of  the  data  for  the  solubility  of  Pb  in  liquid  Zn. 
**&<*».<)  *  12,900  ±  1300;  dS =  7.0  ±  2.0. 


Lead-Copper 

The  solubility  of  Pb  in  solid  Cu  is  low  and  is  neglected  in  calculating  values  of 
S'(Cu  in  Pb.-t).  Excellent  data  are  available  for  a  wide  range  of  temperatures  and  extend 
to  solubilities  as  high  as  NCu  =  0.25.  The  data  of  Kleppa  and  Weil  (8),  Pelzel  (32),  and 
Heycock  and  Neville  (33)  axe  plotted  in  Fig.  12.  An  excellent  straight  line  is  defined  by 
the  data  of  Kleppa  and  Weil,  and  of  Heycock  and  Neville;  this  line  corresponds  to 

=  6900  *  400 
I'd)  1-6  ±  0.4  . 
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Prepare  2  copies  of  this  route  sheet  and 
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ORIG.  IDENT.  SYMBOL  (Mail  Room  Fill  in} 


•PURPOSES 


1.  FOR  INFORMATION 

2.  FOR  APPROVAL 

3.  PREPARE  REPLY 

4.  PREPARE  ENDORSEMENT 

5.  FOR  NECESSARY  ACTION 

6.  FOR  SIGNATURE 

7.  RETAIN  ENCLOSURES 

8.  RETAIN  COPY 
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